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The broad band spectral properties of galactic X-ray binary pulsars 
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BeppoSAX observed several galactic binary X-ray pulsars during the Science Verification Phase and in the 
first year of the regular program. The complex emission spectra of these sources are an ideal target for the 
BeppoSAX instrumentation, that can measure the emission spectra in an unprecedented broad energy band. 
Using this capability of BeppoSAX a detailed observational work can be done on the galactic X-ray pulsars. In 
particular the 0.1-200 keV energy band allows the shape of the continuum emission to be tightly constrained. 
A better determination of the underlying continuum allows an easier detection of features superimposed onto it, 
both at low energy (Fe K and L, Ne lines) and at high energies (cyclotron features). 

We report on the spectral properties of a sample of X-ray pulsars observed with BeppoSAX comparing the 
obtained results. 

Some ideas of common properties are also discussed and compared with our present understanding of the emission 
mechanisms and processes. 



1. Introduction 

The instrumentation aboard BeppoSAX 
is particularly well suited to study the X-ray 
emission from X-ray pulsars. This class of sources 
is composed by binary systems in which a magne- 
tized rotating neutron star accretes matter from 
a less evolved mass-donor star. The mass-donor 
may be a OB supergiant star as in the case of Vela 
X-l, a Be main-sequence or near-main-sequence 
star as in the case of transient recurrent pulsars 
like A0535+26, a low mass star as in the case of 
4U1626-67. The type of mass donor star strongly 
affects the temporal behaviour on the medium 
(days) to long (years) time scales. The tran- 
sient behaviour is almost completely restricted to 
the subclass of X-ray pulsars that have Oe or Be 
counterparts. 

BeppoSAX has observed some persistent pul- 
sars and one transient pulsar during the first year 
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of its operative life. We report results from the 
observations of some of these sources, emphasiz- 
ing the commonalities and the differences. In par- 
ticular we discuss the observational evidence on 
cyclotron line feature, comparing the observed re- 
sults, also in terms of possible correlations, with 
the expected ones on the basis of the available 
theoretical models. 



2. Observations 

X-ray pulsars are a relevant section of the 
BeppoSAX Core Program || devoted to compact 
galactic sources. During the Science Verification 
Program (SVP) and the AOl regular program a 
number of sources in this class were observed. 
In this paper we report the results on Her X-l 
(SVP), Vela X-l (SVP), 4U1626-67 (SVP), Cen 
X-3 (AOl) and GS1843+00 (AOl). The SVP ob- 
servaton of Cen X-3 will not be discussed here. 
A log of the BeppoSAX observations is reported 
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Table 1 



BeppoSAX log of observations 



Source Name 


Observation Date 


Live Time on Source 






(seconds) 


Her X-l 


1996/07/24-28 (SVP) 


90000 


Vela X-l 


1996/07/14 (SVP) 


21600 


4U1626-67 


1996/08/06 09-10 (SVP) 


97000 


Cen X-3 


1997/02/27-28 (AOl) 


20000 


GS1843+00 


1997/04/04 (AOl) 


22000 



in Table 1 

Details of the single observations can be found 
elsewhere (Her X-l: J§; Vela X-l: jl; 4 U1626- 
67: §; Cen X-3 0,y|; GS 1 843+00 

All the observations we report were performed 
with all BeppoSAX telescopes, covering the en- 
ergy band 0.1-200 keV. The observations were 
performed before the failure of the MECS 1 HV 
module occurred in May 1997, therefore all the 
three MECS telescopes are available. The live 
time in column three is calculated for the MECS 
telescopes. The LECS live time is ~ 40% of that 
due to the fact that the instrument is operated 
only during satellite night time. The HPGSPC 
and PDS live time is ~ 50% of that, due to the 
background measure with the rocking collimators. 

In all the observations but for Vela X-l the 
rocking instruments (HPGSPC and PDS) were 
using the default collimator laws: one rocking ev- 
ery 96s. In the Vela X-l observation the dwell 
time was 50s. All data were transmitted using the 
default direct telemetry modes, apart from the 
observation of Her X-l, in which the HPGSPC 
data were transmitted using the full diagnostic 
mode. 

Here we report and compare the results on the 
pulse-phase averaged spectra in 0.1-200 keV. In 
four of these five sources, we have evidence of the 
presence of a cyclotron line. 



3. Results 

The modeling of the broad-band spectra of 
X-ray pulsars is quite difficult. The Her X-l 
spectrum is a prototype of this complexity [J. 
It needs several different spectral components to 



be reasonably modeled: a low energy cutoff, a 
low energy black body - likely due to reprocess- 
ing of primary X-ray radiation at the magnc- 
tospheric boundary-, a Fe L line, a Fe K line, 
a power law, a high energy exponential cutoff 
and a cyclotron feature. A formal fit to the ob- 
served data has 15-19 free parameters. Only the 
BeppoSAX capability to observe the energy spec- 
trum in more than four decades in energy allows 
the clear separation of the different components 
resulting in formal fits in which the correlations 
between parameters are reasonably low. 

Currently there is no theoretical model that 
provides a parametrized model function to fit the 
observed spectra. Therefore the fit to the ob- 
served count rate spectra can be performed only 
using ad hoc analytical models. The first and 
more widely used model function to describe the 
broad band continuum is the power-law-plus- 
cutoff [^U that gives a quite reasonable empirical 
description of the shape of the energy spectra of 
X-ray pulsars. 

E~ a E < E cuto ff 

E-«xeM- E E f E T° ff ) E = E ™toff (1) 

folding 

This model function has a discontinuity at 
E=E cuta ff. 

Other widely used functions are the Fermi- 
Dirac CutOff (FDCO - El) 



E~ a - 



1 



1 + exp( 



E - E cutoff s 



(2) 



f olding 

and the Mihara Negative and Positive power-law 
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Table 2 

Cyclotron lines in the spectra of X-ray pulsars 



Source Name 


ECycl Feature 


FWHM Cycl. 




keV^ 


keV ( Q ) 


Cen X-3 


28.5 ± 0.5 


6.3±2.0 


4U1626-67 


36.5±1.0 


7±2.8 W 


Her X-l 


42.1±0.3 


14.7±0.9 


Vela X-l 


57.9±1.0 


24.0±1. 


A0535+26 




5(W 



(a) Using a multiplicative Gaussian absorption function 

(b) The FWHM is significantly larger using a Lorenzian function 

(c) Using a Lorenzian function |16[ 



plus Exponential (NPEX - |15|) 



(AiE~ ai + A 2 E +a2 ) x exp(-jj) (3) 

We tried all the model functions listed above for 
the sources in Table |[ All these model functions 
show the same practical problem: they do not 
adequately describe the shape and the bending of 
the cutoff, especially in the energy interval where 
the power law joins the exponential tail. 

We also tried for all sources a "modified" 
power-law-plus-cutoff model, in which we used 
a broken power law: 



E -a 2 
E -a 2 



x exp( 



E-E, 



cutoff ■ 



E < Ej 
E = Ej 
E = E, 



break 



break 



(4) 



cutoff 



E folding 

This empirical function seems to provide the 
best performance in fitting the count rate spec- 
tra on the average, therefore it can be used to 
compare the results from the different sources. 

In Table | we summarize the results of the 
broad band fits concerning the cyclotron line. 
The values reported in this table come from 
fits using the model functions in equation [| to 
describe the broad band continuum. The re- 
sults for A0535+26 are taken from Grove et al. 
16 1 . In this table all measurements apart that of 
A0535+26 were done with BeppoSAX . 

In GS1843+00, the source in our sample with 
the hardest spectrum [jl2|], we do not detect any 
cyclotron feature in the pulse-phase averaged 
spectrum. 



The results in this table suggest that there 
may be a correlation between the cyclotron fea- 
ture centroid energy and its FWHM. This cor- 
relation is reported in Figure 1, where we also 
added the data on the cyclotron feature observed 
in A0535+26 by OSSE @. We stress that the 
FWHM parameter in the fit to the count rate 
spectra shows a possible dependence on the choice 
of the function modeling the broad band con- 
tinuum, therefore this correlation must be taken 
with caution. We tried to limit this effect due 
to the choice of the continuum model function, 
using the same broad band continuum in all fits. 
The data follow a general linear trend (with slope 
a ~ 1.7) with some evident scatter. 

We also caution that the line energy reported 
here may not be the energy of the fundamental 
cyclotron feature, but its first harmonic. This 
may be the case of Vela X-l and A0535+26. 

A summary of the cyclotron lines measured by 
BeppoSAX is shown in Figure 2. In this Figure 
we plot a "modified Crab ratio" as described in 
fUf]. To obtain this ratio we 

• divided the count rate spectrum by the 
measured Crab count rate spectrum 

• multiplied the result by the slope of the 
Crab Nebula photon spectrum (E~ 2,1 ) 

• divided the result by the continuum model 
function used for the broad band fit 

This procedure enhances the deviations from the 
broad band continuum and gives a way to visually 
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compare the different features observed in the dif- 
ferent sources. 

4. Discussion 

The correlation showed in Figure 1 was "qual- 
itatively" predicted by Meszaros and Nagel pcfl 
(see also |l7| ) . The model predicts a width of the 
cyclotron feature proportional to its energy and 
to the square root of the electron temperature of 
the atmosphere 

Acjb — wb I 8 X ln(2) x ^ |cos0| (5) 

In this equation Alob is the line width, is the 
cyclotron line frequency, T e is the electron tem- 
perature and 9 is the viewing angle with respect 
to the magnetic field axis. A better insight on the 
properties of the cyclotron lines can be obtained 
with pulse-phase resolved spectroscopy, as equa- 
tion H suggests that there may be a dependence 
on the viewing angle of the observed line width 
(see also Q). 

However Araya and Harding (1996) (l^] cau- 
tion that, in the limit of a single scattering, the 
line width is not related to the electron tempera- 
ture. 

This ambiguity in the interpretation of these 
observational data points out the need of a more 
detailed and quantitative model for the line prop- 
erties and for the broad band continuum emission 
of X-ray pulsars 
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Figure 1 . Cyclotron feature F WHM versus centroid energy. The dashed line shows a simple least square 
linear fit. Even with some scatter present in data (see text for a discussion), a correlation between these 
quantities seems present. 
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Figure 2. Ratios between the count rate spectrum of the different sources and the Crab count rate 
spectrum. The ratios are normalized to the shape of the broad band continuum (see text for details). 
The shape and position of the cyclotron features are well evident. The spectrum of GS1843+00 does not 
show any deviation from the shape of the broad band continuum 



